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ABSTRACT 

Supernova remnants (SNRs) are a complex class of sources, and their heterogeneous nature has 
hindered the characterization of their general observational properties. To overcome this challenge, 
in this paper, we use statistical tools to analyze the Chandra X-ray images of Galactic and Large 
Magellanic Cloud SNRs. We apply two techniques, a power-ratio method (a multipole expansion) and 
wavelet-transform analysis, to measure the global and local morphological properties of the X-ray line 
and thermal emission in twenty-four SNRs. We find that Type la SNRs have statistically more 
spherical and mirror symmetric thermal X-ray emission than core-collapse (CC) SNRs. The ability 
to type SNRs based on thermal emission morphology alone enables, for the first time, the typing 
of SNRs with weak X-ray lines and those with low resolution spectra. Based on our analyses, we 
identify one source (SNR G344.7— 0.1) as originating from a CC explosion that was of unknown origin 
previously; we also confirm the tentative Type la classifications of G337.2— 0.7 and G272.2— 3.2. 
Although the global morphology is indicative of the explosion type, the relative morphology of the 
X-ray line emission within SNRs is not: all sources in our sample have well-mixed ejecta, irrespective 
of stellar origin. In particular, we find that 90% of the bright metal-line emitting substructures are 
spatially coincident and have similar scales, even if the metals arise from different burning processes. 
Moreover, the overall X-ray line morphologies within each SNR are the same, with <6% differences. 
These findings reinforce observationally that hydrodynamical instabilities can efficiently mix ejecta in 
Type la and CC SNRs. The only exception is W49B, which can be attributed to its jet-driven/bipolar 
explosive origin. Based on comparative analyses across our sample, we describe several observational 
constraints that can be used to test hydrodynamical models of SNR evolution; notably, the filling 
factor of X-ray emission decreases with SNR age. 

Subject headings: methods: data analysis — supernova remnants — techniques: image processing — 
X-rays: ISM 



1. INTRODUCTION 

Supernova remnants (SNRs) are a diverse class of ob- 
jects that play an essential role in the Universe, in- 
cluding driving the dynamics of the interstellar medium 
(IS M) and producin g and distributing most of the met- 
als (jFukugita fc Pee bles 2004) . The morphology and dy- 
namics of young SNRs depend on the distribution of 
the ambient medium and on the structure of the stel- 
lar ej ecta. Self-simila r, spherically-symmetric solutions 

[ exist ()ChevalieHll982|) . and they are used widely to in- 
terpret observational data of young SNRs. However, the 

' ejecta are subject to hydrodynamical instabilities that 
preclude a self-similar description of the expansion, and 
thus the use of hydrodynamical models is necessary. 

A major difficulty at present is bridging these hydrody- 
namical models with observations of young SNRs. A few 
direct observables of individual sources (e.g., expansion 
rates) can be compared easily to theoretical predictions. 
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However, the complexity and heterogeneous nature of 
SNRs has limited the ability to define the observed prop- 
erties of SNRs as a class. As a consequence, previous 
observational SNR work has largely focused on interpre- 
tation of single objects, without systematic comparison 
between sources (with some exceptions: e.g., Badenes 
et al. 2010; Long et al. 2010). Although each SNR is 
unique and complicated when studied in detail, it is vital 
to unify the observed characteristics of SNRs to test and 
to improve hydrodynamical models of their dynamics. 

With the advent of high-resolution, space-based tele- 
scopes in the last couple decades, the time is ripe to 
undertake this task. In particular, the Chandra X-ray 
Observatory has facilitated an unprecedented view of 
young ejecta-dominated SNRs since its launch in 1999. 
The sub-arcsecond spatial resolution and the spatially- 
resolved spectroscopy capabilities of Chandra have facil- 
itated detailed studies of the metal-rich ejecta from SN 
explosions as well as their interactions with the surround- 
ing as they expand (see reviews by Weisskopf & Hughes 
2006; Badenes 2010). Chandra has observed over one- 
hundred SNRs in the Milky Way galaxy fCrecn 200^ 
and many others in nearby galaxies (e.g., M33: Long et 
al. 2010). This wealth of data provides the necessary ba- 
sis to characterize the observed X-ray properties of SNRs 
as a class. 

Toward this end, in this paper, we use quantitative 
methods to examine the Chandra images of all SNRs 
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TABLE 1 
Sources, Sorted by Age 



Number Source ObsID ACIS Exp. Age'' Distance Radius'^ Lx'^ References 

(ks) (years) (kpc) (pc) (xlO'^'' erg s~^) 



Type la Sources 



1 


0509-67.5 


776, 7635, 8554 


113 


350-450 


50 


5.96 


1.76 


1 


2 


Kepler 


6714-6718, 7366 


751 


405 


5.0 


3.88 


0.19 


2 


3 


Tycho 


3887 


150 


437 


2.4 


3.72 


0.12 


3 


4 


0519-69.0 


118 


40 


400-800 


50 


6.56 


1.06 


4 


5 


N103B 


125 


37 


860 


50 


5.96 


1.67 


5 


6 


G337.2-0.7 


2763 


49 


~750-10000 


8.0 


7.63 


0.71 


6 


7 


DEM L71 


775, 3876, 4440 


148 


~4360 


50 


11.9 


0.77 


7 


8 


0548-70.4 


1992 


60 


~7100 


50 


17.9 


0.96 


8 


9 


G272. 2-3.2 


9147, 10572 


65 


~6250-15250 


5.0 


13.12 


0.05 


9 


10 


0534-69.9 


1991 


60 


~10000 


50 


21.47 


0.97 


8 


Core-collapse Sources 


11 


Cas A 


4634-4639, 5196, 5319-5320 


993 


309-347 


3.4 


3.77 


2.58 


10 


12 


W49B 


117 


55 


~1000 


8.0 


6.30 


4.48 


11 


13 


G15. 9+0.2 


5530, 6288, 6289 


30 


~1000 


8.5 


8.72 


1.38 


12 


14 


Gil. 2-0.3 


780-781, 2322, 3909-3912 


95 


1623 


5.0 


4.17 


1.18 


13 


15 


Kes 73 


729 


30 


500-2200 


8.0 


6.11 


0.94 


14 


16 


RCW 103 


970 


49 


~2000 


3.3 


5.43 


2.21 


15 


17 


N132D 


5532, 7259, 7266 


90 


~3150 


50 


21.47 


9.92 


16 


18 


G292.0+1.8 


6677-6680, 8221, 8447 


516 


~3300 


6.0 


9.45 


0.56 


17 


19 


0506-68.0 


2762 


38 


~4600 


50 


17.89 


2.62 


18 


20 


Kes 79 


1982 


30 


~6400 


7.1 


13.55 


0.10 


19 


21 


N49B 


1041 


35 


10000 


50 


23.85 


2.39 


20 


22 


B0453-685 


1990 


40 


13000 


50 


20.27 


0.54 


21 


23 


N206 


3848, 4421 


69 


~25000 


50 


29.82 


0.13 


22 


Unknown Type 


24 


G344.7-0.1 


4651, 5336 


27 




14.0 


18.37 


4.23 


23 



Ref erences. — [1] IBadenes etaLl ([20081); [2] [Reynolds eTaLl l [2ro7| ) :J31 1 War ren etHI ll2005t): [41 IRest et al.| l( 2005f): [5] ILewis et aLl 
ll2003) ; [6] [Rakowski et al.' (2006); [7] Hughes et al. (2003); [8] Hendrrck et al. (2003); [9] Harrus et al. (2001); [10] Hwang et aTl 
(pool) : [11] i Lopez et a l. (2009b); [12] Reynolds et al. (2006); [13] Kaspi et al. (2001); [14] Gotthelf fc Vasisht' ( 1997) ; [15 ] Carter etJaP 
aMf): [161 TBorkowski et al.l ( 1200^1); [171 iPark et al.l (120071): [ 181 Hughes et ah (120061 '); [19] iSun et ahl (|200|); [20] fPark et al.l mWI )- 
[21] G aensler et al.l (|2003l ): pb] IWilliams et al.l (120051 ^ [23] lYaina uchi et al. ( 200^ 
^ Aside from historical SNRs and objects with detected light echoes, these ages are very uncertain. 

Radius R, selected to enclose the entire source in the full-band X-ray image, and it is determined assuming the distances listed above. 
" X-ray luminosity in the 0.3-2.1 keV band, from the Chandra SNR Catalog. The v alues for G15.9+0.2, G272. 2-3.2, and G344.7-0.1 is 
calculated using an absorbed planar shock (vpshock) model of the integrated spectra (IRevnolds et al.|l2006l) . 



with strong X-ray line and thermal emission in the Milky 
Way and Large Magellanic Cloud (LMC). We apply well- 
established mathematical tools to characterize the global 
and local morphological properties of SNRs; the statisti- 
cal approach we take here enables, for the first time, the 
capability to compare our results within and between 
sources to set observational constraints on hydrodynam- 
ical models. 

This paper is organized as follows. In §2, we describe 
the sample and observations used in this study, and in 
§3, we present the methods employed in our analyses. 
In §4, we give our results, and in §5, we examine how 
morphological properties vary across our sample. §6 dis- 
cusses the observational constraints on hydrodynamical 
models we have found in this paper and the implications 
of our findings. 

2. OBSERVATIONS AND DATA PREPARATION 

For our analyses, we utilize archival Chandra ACIS ob- 
servations of the twenty-four SNRs listed in Table 1 and 
shown in Figure [1] Ten of our sources are thought to be 
from Type la SNe, thirteen are considered to have origi- 
nated from CC SNe, and one is unknown (see Table 1). 



We selected sources in the Milky Way and LMC that 
have prominent thermal emission from ejecta in their 
global X-ray spectra (with 0.5-2.1 kcV X-ray counts per 
unit area >0.01 counts/pixel^ within the radius R that 
encloses their signal in Table 1). This criterion excludes 
SNRs dominated by n on-thermal emission, like Gl. 9-1-0. 3 
(jRevnolds et al.|[2008() . We also required that the sources 
be imaged fully in one ACIS pointing; this restriction re- 
moved SNRs with large spatial extent, such as SN 1006 
(e.g., Long et al. 2003). Additionally, we exclude SNRs 
that are interacting with or are distorted by molecular 
clouds (e.g., N63A: Chu & Kennicutt 1988, Warren et 
al. 2003; G349.7-h0.2: Lazendic et al. 2005) and those 
whose large-scale morphologies are the result of axisym- 
metric winds from pulsars (such as 3C 58, Slane et al. 
2004). 

Each source was observed for ~30-1000 ks. Data re- 
duction and analysis was performed using the Chandra 
Interactive Analysis of Observations (ciAO) Version 4.0. 
We followed the CiAO data preparation thread to repro- 
cess the Level 2 X-ray data, and we extracted global 
X-ray spectra of each source using the CIAO command 
specextract. Then, we produced exposure-corrected im- 
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Fic;. 1. Chandra X-ray soft-band (0.5-2.1 kcV) images of the 24 SNRs listed in Table 1. The cyan circles mark the full-band centroids 
of each SNR used in our power-ratio/multipole expansion method. Numbers correspond to those in Column 1 of Table 1. Red numbers 
denote Type la SNRs; light blue numbers denote CC SNRs. 



ages of the soft X-rays (0.5-2.1 keV; we set the upper 
limit to 2.1 keV in order to include all of the flux in 
the Si XIII line) and emission lines (listed in Table 2) 
for each source. For sources with bright pulsars (like 
Gll.2-0.3, Kes 73, RCW 103, and G292.0-H1.8), the 
pulsar location and extent was identified using the CIAO 
command wavdetect (a source detection algorithm using 
wavelet analysis; Freeman et al. 2002). We replaced 
the region identified by wavdetect with pixel count val- 
ues selected from the Poisson distribution of the area 
surrounding the pulsar using the CiAO command dmfilth. 
This process removed the bright pulsars while preserving 
the morphologies of the diffuse emission surrounding the 
pulsars. Generally, the removed area of the pulsars was 
small (^ 16 pixels^). In sources where the pulsar emis- 
sion was more extended (e.g., RCW 103 and Kes 73), 
it was necessary to replace a larger area (^ 400 pixels^); 
since these SNRs are shell-like (Green 2009), this pro- 
cedure did not alter their overall morphology. No other 
point sources were removed because of potential confu- 
sion with small ejecta substructures. 

Given the young-to-middle age of our sources (see Ta- 
ble 1), we expect that all are ejecta-dominated, and the 
shocked ISM has only a minor contrib ution to the ob- 
served X-ray flux (jBadenes et al.l[2010l ). X-ray spectral 
modeling generally confirms that abundances are above 
those of the ISM (see references in Table 1), suggest- 
ing the emission is indeed dominated by the shocked SN 



ejecta and not by shocked ISM. 

3. METHODS 

In what follows, we apply two methods to quantify 
the X-ray morphologies of our SNR targets: a power- 
ratio method (PRM) to measure symmetry and wavelet- 
transform analysis (WTA) to probe X-ray substructure. 
These techniques were introduced in Lopez et al. 2009a 
(L09a, hereafter); we refer the reader to that paper for a 
detailed formalism. Here, we give a brief overview of the 
methods and their uses. 

3.1. Power-Ratio Method 

The PRM enables the measurement of asymmetries in 
X-ray surface brightness distributions, and we employ 
this technique here to compare the global morphologies 
of thermal emission in Type la and CC SNRs. The 
method was first applied to characterize the X-ray mor- 
phology of galaxy clusters observed with ROSAT (Buote 
fc Ts ai, 1995, 1996) and with Chandra (jJeltema et al.l 
120051 ). Subsequently, L09a and Lopez et al. 2009b (here- 
after, L09b) developed and extended the technique to 
Chandra observations of SNRs. The PRM measures 
asymmetries in an image via calculation of the multi- 
pole moments of the X-ray surface brightness in a cir- 
cular aperture. It is derived similarly to the multipole 
expansion of the two-dimensional gravitational potential 
within an enclosed radius R: 
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(j)) = -2Gaa In ( | - 2G 



— i^ra COS mcj) + bra sin mcj)) , (1) 



where the moments am and are 

am{R)^ I Y.{x'){R'y" COB met)' (fx\ 

J R'<R 

iim{R) = / (i?')™ sinm^'d^x', 

x' — {R' , (p'), and S is the surface mass density. For our 
imaging analyses, the X-ray surface brightness replaces 
surface mass density in the power ratio calculation. 

The powers of the multipole expansion are obtained 
by integrating the magnitude of ^'m (the mth term in 
the multipole expansion of the potential) over a circle of 
radius R, 

1 f^^ 

^r„(i?) = ^y (2) 
Ignoring the factor of 2G, this equation reduces to 



Po-[ao In {R)Y 
1 



(3) 



The moments am and bm (and consequently, the pow- 
ers Pm) are sensitive to the morphology of the X-ray 
surface brightness distribution, and higher-order terms 
measure asymmetries at successively smaller scales rela- 
tive to the position of the aperture center (the origin). 
To normalize with respect to flux, we divide the powers 
by Pq to form the power ratios, Pm/Po- Pi approaches 
zero when the origin is placed at the surface-brightness 
centroid of an image, so we have set the aperture center 
in all analyses to the full-band (0.5-8.0 keV) centroid of 
each remnant. In this case, morphological information is 
given by the higher-order terms. P2/ Po is the quadrupole 
ratio; examples of sources that have high P2 / Po are those 
with elliptical/elongated morphologies or those with off- 
center centroids because one side is substantially brighter 
than the other. P3//0 is the octupole ratio; examples of 
sources that have high P3/P0 are those that have asym- 
metric or non-uniform surface-brightness distributions. 

A Monte Carlo approach described in L09b is used to 
estimate the uncertainty in the power ratios. Specifi- 
cally, the exposure-corrected images (normalized to have 
units of count s) are adaptively-binned using the program 
AdaptiveBin (jSanders fc Fabianll2001|) such that all zero 
pixels are removed, smoothing out noise. Then, noise 
was put back in by taking each pixel intensity as the 
mean of a Poisson distribution and selecting randomly a 
new intensity from that distribution. This process was 
repeated 100 times for each soft-band image, creating 
100 mock images per source. The 1-a confidence limits 
represent the sixteenth highest and lowest power ratio 
obtained from the 100 mock images of each source. 



3.2. Wavelet- Transform Analysis 

Wavelet-transform analysis (WTA) is the other 
method we use in this paper to characterize the X-ray 
morphology of SNRs. In L09a, we demonstrated that 
this technique can measure accurately the substructure 
and filling factor of X-ray emitting plasma, and we ap- 
plied the method to one complex SNR, W49B. Below, we 
employ WTA to compare the X-ray substructure prop- 
erties of the line emission in our SNR sample. 

WTA was first applied successfuly to RO S AT a.nd Ein- 
stein data to extr act the small-scale X -ray structure 
of galaxy clusters ()Grebenev et al.l [19951 ). A wavelet- 
transformed image is a decomposed image of a signal's 
intensity (from herein, power) measured at the scale of 
a filter size. Mathematically, a wavelet transform w is 
the correlation of a signal s(x, y) in an image with the 
analyzing wavelet function g(x, y): 



I \ I \ ^ ( ^ y 

w(x, y, a) = s{x, y) -5 - 
a \a a 



(4) 



where a is the scale (or width) of the wavelet transform. 
We utilize a radial Mexican-hat function g{-,-) (the nor- 
malized second derivative of a Gaussian function) of the 
form 

g^^^y^^(2^^^\^-i.^W)/2a^, (5) 

a a \ a'' J 

Wavelet-transformed images are produced by calculat- 
ing w^x, y, a) for each pixel (m, n) in a raw image: 



w{m, n, a) 



-2^c,jg[ , ), (6) 



where is the number of counts in the pixel. 

Essentially, w measures the summed intensity enclosed 
by the area of the Mexican hat. Thus, the size of an 
individual source can be characterized by the scale where 
the convolution of the wavelet and a signal reaches a 
maximum. The wavelet transformation of an isotropic 
Gaussian signal of size a and intensity / at a position 
(xo.yo) is 



w{xo,yo,a) 



21 



1 



(7) 



If we divide this relation by a, it has an absolute max- 
imum at a — a, which we define as amax- Thus, by 
identifying the peak in a plot of w/a versus a for the 
central pixel of an emitting substructure, we can mea- 
sure its size. In addition to measuring the scale of in- 
dividual substructures, we can sum all the pixels in the 
wavelet-transformed images at each scale a to find the 
power profile of an entire source. The resulting power 
profile, (w) I a versus a, depends on the scale of isolated 
structures as well as the filling factor of the emitting ma- 
terial. 

To aid in understanding the method, we describe here 
how the power profiles would look for a variety of cases. 
These examples are given quantitatively in §2.3.2 of 
L09a. In an image with noise only, the power profile 
would peak at scales of a single pixel and decline rapidly 
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Fig. 2. — Left: Power ratios, the quadrupole ratio P2/P0 versus the octupole ratio P3/P0, of the soft X-ray band (0.5-2.1 keV) for 
twenty-four SNRs in the Milky Way and LMC. Right: The same plot using only Si XIII 75-2.0 keV) in the seventeen SNRs from 
L09b. Type la SNRs are in red, CC SNRs are in blue, and 0548—70.4 is in purple because of its anomalous ejecta properties that make its 
type uncertain. The quadrupole ratio is a measure of ellipticity/elongation, and the octupole ratio quantifies the mirror asymmetry of the 
emission. We find that the Type la SNRs are more circular and symmetrical than the CC SNRs. 



toward zero. If an image has only one substructure with- 
out noise, the power profile would have a global maxi- 
mum at the scale of that substructure. With noise, the 
power profile should be identical to that of the no-noise 
case, as long as the signal-to-noise ratio is greater than 2 
and the noise-only pixels (those with Cmax = 1 pixel) are 
removed. As the number of substructures in an image in- 
creases, their emission will agglomerate, augmenting the 
surface area of the emitting regions (i.e., increasing the 
filling factor) and causing Omax of the power profile to 
increase. 

4. RESULTS 

We use the methods from §3 to the sample in Table 
1 to examine the global and local X-ray morphological 
properties of SNRs. 

4.1. Global X-ray Morphologies 

To measure the global X-ray morphologies of the ther- 
mal emission in SNRs, we applied the PRM to the soft X- 
ray (0.5-2.1 keV) images of the twenty-four SNRs shown 
in Figure [T] This work is an extension of the analyses 
in L09b, where we employed the PRM and calculated 
the multipole moments of the Si xiii images of seventeen 
Galactic and LMC SNRs observed by Chandra. In L09b, 
we found that the CC and Type la SNRs can be distin- 
guished by their quadrupole and octupole ratios, P2/P0 
and P3/P0 respectively. In particular, the CC SNRs had 
an order of magnitude greater P2 / Pq than the Type la 
SNRs, indicating CC SNRs are statistically more elon- 
gated/elliptical than Type la SNRs. Additionally, the 
CC SNRs had a factor of two larger P3/P0 than the Type 
la SNRs, suggesting the CC SNRs are more mirror asym- 
metric than Type la SNRs. The results were the same 
for other X-ray emission lines besides Si Xiii, e.g. Ne ix, 
Mg XI, and S xv. 



Here, we apply the method to thermal X-rays in SNRs 
generally. In doing so, we are able to increase the sample 
size since several remnants have strong bremsstrahlung 
emission without resolved or strong emission lines. Thus, 
in addition the seventeen targets from L09b, seven 
new sources have sufficient bremsstrahlung emission 
for our analyses: G337.2-0.7, G272.2-3.2, 0534-69.9, 
0506-68.0, Kes 79, N206, and G344.7-0.1. To ensure 
that we are measuring thermal X-rays and not non- 
thermal emission, we analyzed the soft X-ray images de- 
scribed in §2, since bremsstrahlung dominates over syn- 
chrotron emission below ^2 keV. 

Figure 1 (left) shows the resulting P2/P0 versus P3/P0 
plot for the soft X-ray images; the analogous Si XIII plot 
from L09b is given (right) for comparison. We find that 
the CC SNRs have a mean P2/P0 = (94.2±0.4) x lO"'^ 
with a standard deviation of 90.4xl0~^, and the Type 
la SNRs have a mean P2/P0 = (6.53±0.05) x lO"'^ with 
a standard deviation of 7.89 xlO~^ (excluding SNR 
0548-70.4, see the discussion in L09b). The mean P3/P0 
of the two classes are also different: the mean of the Type 
la SNRs is (2.60±0.13)xl0-'^ with a standard deviation 
of 1.90 and of the CC SNRs is (5.01±0.88)xl0~^ with a 
standard deviation of 4.33. This discrepancy in P3/P0 
can be attributed to the CC SNRs with large P3/P0 
(J; 60). Generally, our findings are consistent with those 
of L09b: CC SNRs are much more asymmetric or ellip- 
tical than Type la SNRs. We attribute these differences 
to the distinct explosion mechanisms and circumstellar 
medium structures of Type la and CC SNRs. 

Of the Type la SNRs, Kepler has one of the largest 
P2 1 Pq because of its off-center centroid (see Figure [T]) 
since one side being brighter than the other. Sources 
with more symmetric and homogeneous emission (like 
G272.2-3.2) have the smallest P2/^b- Additionally, 
centrally- filled SNRs (e.g., N103B) tend to have smaller 
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Fig. 3. — Chan dra ACIS full-band (0.5—8.0 keV) spectra for the nine sources analyzed in §4.2. Black lines indicate which emission features 
we analyzed in § 14.21 these emission lines are listed in Table 2, from low to high energy. 



P2/PQ as well. Of the CC SNRs, the sources with 
bright pulsars tend to have the lowest P2/P0 (such as 
B0453— 685 and Kes 79) suggesting those SNRs are more 
circular and symmetric than those without pulsars or 
neutron stars. Finally, SNRs with elongated or elliptical 
shapes (like W49B) have the highest P2/P0, and those 
with large-scale asymmetries have the greatest P3/P0 
(e.g., ROW 103). 

4.2. Small-Scale Structure 

From § 14.11 it is evident that the large-scale morpholog- 
ical differences of the X-ray line and the thermal emitting 
material between SNRs can be used to distinguish the 
explosion type. Next, we consider the relative morpholo- 
gies of different X-ray lines within individual sources and 
what their properties can reveal about their explosions 
and dynamical evolution. The results of §4.1 hint that 
the X-ray lines of each SNR have similar morphologies, 
since the PRM could predict accurately the explosion 
type regardless of which line image was analyzed. 

Toward this end, we apply the WTA technique out- 
lined in §3.2 to all the emission line images of our sources. 
Since we are focusing on the comparison between emis- 
sion lines, we limited our sample to only those with at 
least two strong X-ray line features (with counts per unit 
area >0.01 counts/pixeP). Additionally, as we are con- 



TABLE 2 

X-RAY Emission Line Selection 



Source 



Lines^ 



Cas A O Cont'', Mg xi, Si xiii, S xv, Ar xvii, Ca xix, Fe xxv 

Kepler Fe L, Mg XI, Si xiii, S xv, Ar xvii, Ca xix, Fe xxv 

Tycho Fe L, Si XIII, S xv 

W49B Si XIII, S XV, Ar xvii, Ca xix, Fe xxv 

G15. 9-1-0.9 Si xio, S xv 

Gil. 2-0.3 Mg XI, Si XIII, S xv 

Kes 73 Mg xi. Si xiii, S xv 

RCW 103 Fe L, Mg xi. Si xiii, 

G292.0+1.8 O viii, Ne ix, Mg xi. Si xiii, S xv 



^ Energy ranges for individual lines vary slightly across the sources, de- 
pending on e.g., the width of the lines. On average, the bands are: O Cont: 
0.6-0.8 keV; O viii: 0.6-0.7 keV; Ne ix: 0.85-0.95 keV; Fe L: 0.9-1.1 keV; 
Mg XI: 1.20-1.50 keV; Si xiii: 1.7-2.1 keV; S xv: 2.25-2.60 keV; Ar xvii: 
2.9-3.3 keV; Ca XIX: 3.7-4.1 keV; Fe xxv: 6.2-6.9 keV. 

The oxygen in Cas A is expected to be completely ionized and to dominate 



the bremsstrahlung continuum jVink et al.|[T996D . Theref 
0.6—0.8 keV continuum as a proxy for the oxygen 



ore, we use 



the 



sidering local structures with arcsecond extents, we re- 
strict the analyses of this section to Milky Way SNRs 
to ensure that we can resolve sub-pc scale structures. 
Thus, we limit our sample to the nine SNRs that satisfy 
these criteria. Figure |3] gives the global X-ray spectra 
for these nine SNRs, with the black lines labeling all the 
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X-ray emission lines whose images we analyzed. Table 2 
lists these X-ray emission lines, in order from the lowest 
to highest energies; the nine SNRs had 2-7 lines that had 
sufficient surface brightnesses for our analyses. 

Figures IMTSl (given in the Appendix §A.l) show the 
raw images of these X-ray lines as well as the result- 
ing wavelet-transformed images at five scales for each 
SNR (except W49B, given in Figure 14 of L09a). Since 
the wavelet essentially acts like a filter to pick up the 
emission at different scales, each field displays the X-ray 
power of the different lines at the given sizes. At small 
scales, noise and random fluctuations dominate, and with 
increasing filter sizes a, the distribution and substructure 
of each ion becomes more evident. Generally, the eight 
SNRs have similar morphologies among all of their emis- 
sion lines. These results contrast the case of W49B: in 
L09a, we found that the iron in W49B was largely ab- 
sent in the Western half of that source, whereas the other 
elements were more symmetrically distributed. 

Using the transformed images, we identify and measure 
the scale rc and position of individual, isolated X-ray sub- 
structures of the line emitting material in each SNR. We 
find 15-45 substructures in each source, and they span a 
range of scales (~3"-35"). We only include substructures 
with scales ^ 3" to ensure that point-spread-function ef- 
fects do not influence the results. 

As a probe of chemical mixing, we can compare the 
substructure sizes Tc and locations of different elements 
within each source. Specifically, we identify and measure 
substructures which are spatially coincident (defined as 
those less than ten pixels w 5" apart) in two X-ray line 
images. When possible, we restrict our analyses to ions 
that are the products of different burning processes, and 
thus the relative scale and position of substructures re- 
veals the effectiveness of chemical mixing in the SNRs. 
For six sources, we utilize the two strongest X-ray emis- 
sion lines. Si xiii (a product of oxygen burning) and Mg 
XI (a product of carbon and neon burning). In Tycho 
(which lacks a prominent Mg XI feature), we compare 
the substructures of Si xiii and S xv (both from oxygen 
burning) as well as Si xiii and Fe L (a product of silicon 
burning). For the other SNRs (W49B and G15.9-H00.3), 
we are limited to Si xiii and S xv, since these SNRs have 
only products of oxygen burning (except for Fe XXV in 
W49B, which has a disparate morphology relative to the 
other ions; L09a). 

Figure m plots the resulting substructure sizes rc in our 
nine sources. We find that ^^90% of identified (i.e., the 
brightest) substructures in one X-ray line image have a 
corresponding substructure in the other image, suggest- 
ing the elements are well mixed throughout the SNRs. 
Broadly, the slopes of the plots in Figure |4] are consis- 
tent with unity, indicating that the substructures of the 
X-ray line emitting material within each SNR have sim- 
ilar physical scales, ranging from 1-16% of the radius of 
each SNR. 

4.3. Average Power over Many Scales 

In addition to comparing individual structures on small 
scales, we can use WTA to examine the overall power pro- 
files (power as a function of scale over the entire source) 
of the ions. Figure Ogives the power profiles {{w)/a 
versus a) for the emission lines in each remnant. The 
curves for each SNR reach maxima at different a„iax, 



ranging from ^10"-50". Since the individual substruc- 
tures identified above are generally smaller than Omax in 
each source, the scale a„iax is a reflection of the surface 
filling factor: the greater the value of flmax, the larger 
the filling factor of the emitting material. 

Within each source, the power profiles of the emission 
lines have similar shape and identical maxima, with only 
a few exceptions. The Fe xxv in W49B peaks at 25% 
larger scales than the other ions (22.5" versus 17.5") of 
that source; this discrepancy is the largest among our 
nine sources. One exception is the O in Gas A, which 
peaks at amax — 25"while the other ions have maxima at 
timax = 30". Another exception is the Si xiii in Kes 73, 
which peaks at slightly larger scales, Omax = 13", than 
the Mg XI and S xv, with a,„ax = 12". Additionally, 
some remnants, such as Gas A, Kepler, and G292.0-t-1.8, 
probably have excess power at small scales in the higher- 
energy lines (like S xv, Ar xvii, and Ga xix) because of 
contamination from non-thermal emission. 

To test whether these differences in the ion power 
profiles are significant, we compare the curves for each 
SNR quantitatively by measuring their cumulative power 
across a range of sizes. In particular, we determine the 
fraction of total power that each ion image has at scales 
above and below their a,nax values. For this analysis, in 
sources where the ions have different Omax, we used the 
flmax of Si XIII. We find that all the SNRs except W49B 
have <6% differences between their ions' relative power 
above and below scales of Cmax- By constrast, the Fe xxv 
in W49B has w34% less power than Si xiii below scales 
of its Omax ~ 17.5". We conclude that variations in the 
SNRs' power profiles are minor (excluding W49B), and 
the elements of our sources have similar surface filling 
factors. 

5. SUBSTRUCTURE TRENDS ACROSS SNRS 

Here, we examine how these substructure character- 
istics vary across SNRs. First, we compare the power 
profiles between SNRs to examine whether they depend 
on age. Toward this end, we plot {w)/a of Si xiii in ev- 
ery source versus scale a in physical units (parsecs), as 
shown in Figure |6] (left). The Omax values vary by only 
a factor of ~3, from wO.24 pc for G15.9+0.2 to «0.67 pc 
for W49B. This value does not appear to depend on age: 
for example, Tycho has roughly the same maximum as 
G292.0-f 1.8, flmax ~ 0.57 pc, even though G292.0+1.8 is 
almost a factor of ten older than Tycho. The uncertain 
distances of Kepler and G15.9-I-0.2 may be the reason 
those two sources are outliers from the other SNRs. 

Since the conversion of scale a to physical units de- 
pends on the sometimes largely uncertain distance to our 
sources, in Figure [6] (right), we also plot the power pro- 
files as a function of the dimensionless quantity a/Rx- 
Here, Rx is the radius of the X-ray emission in the full- 
band image. In elongated sources (such as W49B), we 
define the radius Rx as the semimajor axis of the el- 
lipse that encloses the remnant's full-band X-ray surface 
brightness. The SNRs' maxima a^mx/ Rx span a decent 
range, from a,-ae,x/Rx « 0.04 (G15.9+0.2) to Oraax/ Rx ~ 
0.57 (G292.0-hl.8). 

Broadly, the younger SNRs of our sample have larger 
amax/Rx than the older sources, suggesting that the fill- 
ing factor of the emission decreases with age. Some rem- 
nants (e.g., Tycho and G15.9-I-0.2) do not follow this 
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Fig. 4. — Comparison of individual substructure sizes for the two strongest emission lines (Si XIII and either Mg XI or S xv) in each 
remnant. When available, we compared substructure scales of elements from different burning processes; therefore, Tycho also includes 
the comparison of Fe L versus Si XIII in gray. The error bars reflect the uncertainty in the size estimate, and the dashed lines have slope 
of unity. 90% of identified substructures in one X-ray line image had a corresponding substructure in the other X-ray line image of that 
SNR, and the substructure sizes of the X-ray lines within each source are nearly identical. Only substructures with scales >6 pixels were 
included to avoid PSF effects. 



trend, however. The explanation for this anomaly is un- 
certain, but we note that Tycho and G15.9+0.2 have the 
shortest ionization timescales Uct (defined as the prod- 
uct of the electron density Uc with the time since the 
plasma was shocked) of our nine sources: Tycho has a 
mean rict ^ 3 x 10^" s cm~'^ (based on our spectral anal- 
ysis des cribed below) and G1 5.9-f0.2 has Uct ^ Q x 10^" 
s cm-3 (jRevnolds et al.ll2006[ ). 

Next, we investigate the relationship between individ- 
ual substructure size and luminosity. To obtain the sub- 
structures' emitted fluxes (i.e., absorption corrected), we 
extracted and modeled the Chandra X-ray spectra of ev- 
ery substructure identified in the analysis from Figure 01 
Spectra were extracted from all available observations 
with regions of radii corresponding to the scale identi- 



fied by the WTA. For all the sources except W49B, we 
fit these spectra using an absorbed, variable abundance 
plane-parallel shocked plasma model with constant tem- 
perature, phabs X vpshock, in XSPEC Version 12.4.0. 
Previous X-ray analysis of W49B has shown it is in col- 
lisional ionization equilibrium (GIF) and requires two 
plasmas to sufficiently fit its spectra (Miceli et al. 2006; 
L09a). Therefore, for W49B, we use instead two GIF 
components for this source: one cool plasma with fixed 
solar abundances (XSPFG model mekal) and one hotter 
plasma with varying supersolar abundances (vmekal). In 
the case of Gas A, the continuum emission is thought 
to be from completely-ionized oxygen rather than H and 
He (Vink et al. 1996), so we set the abundances of H and 
He to zero. Since three sources (Gil. 2— 0.3, Kes 73, and 
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Fig. 5. — Relative power versus substructure scale {(to/a) vs. a; in arcseconds and parsecs) for all the emission lines in each remnant (see 
Table 2). The scale a where the power of an individual structure reaches its maximum, Omax, gives its size. The power profiles of ions in 
each source have similar shape, with only < 6% fractional differences in the profiles of a given source. The only exception is W49B, where 
the Fe XXV has substantially less power (Ri34%) at small scales than the other ions. 



RCW 103) do not have published Chandra X-ray spectra 
and modeling of their ejecta, we provide a more detailed 
analysis and discussion of these sources in §A.2. From 
the best-fit spectral models of the nine SNRs, we deter- 
mined the emitted flux of each substructure in the Si 
XIII line (over the range 1.75-2.0 keV) and measured the 
luminosity assuming the distances in Table 1. 

Figure [7] shows the resulting plot of the Si Xiii lumi- 
nosity versus substructure size for the nine SNRs. We 
fit the log-log data with a linear polynomial of the form 
log Lgi = b{\og r^) + c; Table 3 lists the best-fit slopes 
b and y-intercepts c and associated errors of these anal- 
yses. All the slopes h are consistent with values of 2-3. 
We would expect that if the substructures are optically 
thin, they would have a power-law index of 3; therefore, 



we attribute the slopes 1.91 < 6 < 2.97 to a low volume 
covering fraction of the individual unresolved substruc- 
tures. In physical terms, the y-intercept c of these fits is 
the extrapolated luminosity of a substructure that is one 
parsec in size. Our results show that this parameter de- 
creases with age, spanning over two orders of magnitude 
from the youngest (Cas A) to oldest (G292.0-I-1.8) of our 
sources. This trend suggests that individual substruc- 
tures become less luminous with time, consistent with 
the result that filling factor decreases with age. 

We obtain several other physical parameters of the in- 
dividual substructures with the best-fit spectral models, 
including their electron temperatures kT and ionization 
timescales Uct. We searched for trends between kT and 
substructure size as well as Uct and substructure size. 
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Fig. 7. — Si XIII luminosity Lsi (1.75-2.0 keV) versus individual 
substructure size rc for all the emitting regions in the nine SNRs. 
Best-fit lines to this log-log plot are given in Table 3. We find 
that substructures in older sources are less luminous than those in 
younger SNRs. 

and we found no clear relation within nor between SNRs. 
Since these figures are essentially scatter plots, we do not 
reproduce them in this paper. 

6. DISCUSSION 

In this paper, we have exploited the wealth of Chan- 
dra ACIS data on galactic and LMC SNRs to examine 
the observed X-ray properties of these sources. We have 
applied statistical tools to every remnant with strong 
line and thermal emission to enable comparison of the 
local and global morphological characteristics of SNRs. 



TABLE 3 
Best-Fit Lines for Figure [T] 
logLsi = b(logrc) -I- c 

Source b c 



Cas A 
Kepler 
Tycho 
W49B 
G15. 9-1-0.9 
Gil. 2-0.3 

Kes 73 
RCW 103 
G292.0-I-1.8 



2.11^ 



1-0.49 
-0.48 

2.97±0.37 
2.39±0.32 

^•*°-0.66 

2.53±1.87 

,+0.53 
-0.52 

2.55±0.53 

r, r,r, + 0.37 

^•^^-0.36 

1.91±0.81 



2.73^ 



35.8«;^ 

34.9±0'.3 
34.6±0.4 

QQ Q+0.8 

■JJ-O-O.g 
34.7±0.4 
34 4+0-4 

33.6t°* 
33.1±0.6 



Ultimately, we aimed to determine constraints on the 
physical processes underlying the dynamical evolution of 
SNRs. 

In §4.1, we demonstrated that the large-scale mor- 
phologies of the X-ray line and thermal emitting material 
are different for Type la versus core-collapse SNRs: the 
Type la SNRs have statistically more spherical and mir- 
ror symmetric X-ray emission than the CC SNRs. The 
ability to distinguish the explosion type based on the 
bremsstrahlung emission morphology alone enables, for 
the first time, the typing of remnants with weak X-ray 
lines and of those with low resolution spectra. It also sug- 
gests that it may be possible to type SNRs using other 
energies where bremsstrahlung dominates. In our analy- 
sis, we have successfully identified the SNR G344.7— 0.1 
as originating from a CC explosion, a source of un- 
known explosion type previously. Additionally, we have 
confir med the tentati ve classifications of G337.2— 0.7 
(Rak owskiet all I2006D and of G272. 2-3.2 ()Park et al.l 
I2009f ) as Type la SNRs, bringing the number of known 
Type la SNRs in the Milky Way to six (the four others 
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being Gl.9+0.3: Reynolds et al. 2009; Tycho; Kepler; 
and SN 1006). 

In §4.2, we investigated the small-scale structures of 
several ions within nine galactic SNRs. We found that 
the emission lines within the SNRs have remarkably sim- 
ilar substructure locations, scales, and power profiles, 
even if the ions are products of different burning pro- 
cesses. This result implies that the metals within the 
remnants (both in the ejecta and shocked CSM) must 
have similar spatial distributions, and as such, the metals 
within SNRs must be globally well-mixed. These findings 
reinforce observationally that hydrodynamical instabili- 
ties efficiently mix ejecta at the scales we can resolve 
here (although the metals may still be dynamically dis- 
tinct; see Badenes et al. 2005). Our analysis shows that 
these results are true for both Type la and CC SNRs, 
indicating that the mixing efficiency is not dependent on 
explosion type or on CSM structure. Therefore, we con- 
clude that the relative, large-scale morphologies between 
the different X-ray emission lines in a source cannot be 
used to distinguish explosion type. 

From our analyses, we do not find evidence of signifi- 
cant ejecta stratification on the scales studied here ('^3"- 
40") in our Type la SNRs, Kepler and Tycho. Chemical 
stratification is observed in Type la SN (e.g., Mazzali et 
al. 2007), and in some Type la SNRs (e.g., Kosenko et 
al. 2010). Previous X-ray studies of Kepler with XMM- 
Newton have shown that Si K and Fe L have similar 
radial profiles in the north, whereas the Si K ext ends to 
large r radii than Fe L in the south (.Cassam-Chena'i et al.l 
|2004|) . Indeed, at the smallest scales in Figure IHl the 
transformed images of Fe L show some substructures in- 
terior to those of Mg XI and Si Xiii in the south. In the 
case of Tycho, prior work has demonstrated that silicon 
and iron are si milarly distributed, with both close to th e 
forward shock (| Warren et al.ll2005l : iBadenes et al.ll2006f ). 
Therefore, our results are consistent with the emerging 
picture for these two sources: some mechanism (such 
as hydrodynamical instabilities) has reduced the initial 
stratification of the ejecta expected right after their ex- 
plosions. 

From our analyses, the only exception to the above 
statements regarding substructure and mixing is W49B. 
In that SNR, we found that the Fe xxv had distinct mor- 
phological and substructure properties from the lower- 
Z ions (L09a). In L09a, we demonstrated that W49B 
likely originates from a jet-driven/bipolar explosion (e.g., 
Ramirez-Ruiz & MacFadyen 2010), based on these mor- 
phological discrepancies and the abundance ratios and 
masses of the different species. Although Cas A also has 
jet-like features, these structures do not contribute signif- 
icantly to the remnant's overall X-ray surface brightness, 
and the power profiles of all the ions in Cas A are iden- 
tical. Thus, we conclude that only large-scale discrepan- 
cies between ions' spatial distributions within a source 
can facilitate identification of the explosion mechanism. 

Our findings reinforce the unique nature of W49B as 
the remnant of a jet-driven explosion. We can esti- 



mate whether this kind of event is expected to have oc- 
curred recently in the Milky Way (MW) galaxy by con- 
sidering the observed rates of SNe in the Local Group. 
Jet-driven/bipolar explosions are associated with Type 
Ib/Ic SNe, a subclass of CC SNe. I n the MW, the 
rate of CC SNe is ^-^2 per century (jTammann et al.l 
IHqI. Of these CC SNe approximately 20% are Type 
Ib/Ic (jSmartt et al.ll2009D : thus, the rate of Type Ib/Ic 
SNe is ~l/250 years. Some subset of these SNe will 
be jet-driven/bipolar explosions, but it is still uncer- 
tain what fraction of Type Ib/Ic are jet-driven. How- 
ever, the estimated rate of hypernovae (HNe, which are 
super-energetic bipolar SNe) per galaxy is ~ 10~^/year 
(jlzzard et al.ll200l : IPodsiadlowski et al.il2004D . If 1-10% 
of bipolar explosions are HNe, the rate of jet-driven SNe 
would be one every 10000-1000 years. Therefore, it is 
reasonable to expect several jet-driven/bipolar SNe in 
the MW, and a few should be observable as SNRs at 
X-ray wavelengths. 

From the comparative analysis of our sources, we have 
set several observational constraints that are useful tools 
to test the validity of theoretical models of SNR dynam- 
ical evolution. Specifically, we have found that: 

• Mixing of shocked ejecta and CSM must be efficient 
on the scales resolved here (^3-35"); 

• Individual emitting substructures within a source 
span a range of sizes, and ^ 90% of the metals' 
brightest substructures in an SNR should be coin- 
cident and have equal surface emission scales; 

• The surface area of the X-ray emission (given by 
our parameter a„iax) does not depend on age (likely 
because it depends on a combination of age and 
density). However, relative to the size of the rem- 
nant, the surface area of the X-ray emission de- 
creases with age; 

• Individual emitting substructures become less lu- 
minous with time; 

• The scale of individual substructures does not 
tightly correlate with the temperature or ionization 
timescale of that substructure. 

The list above is a first step toward a broad observa- 
tional baseline for direct comparison to the predictions 
of hydrodynamical models. As such, the local and global 
morphological properties described here should aid in ad- 
vancing understanding of SNRs, both from an observa- 
tional and theoretical perspective. 
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APPENDIX 
A.l. WAVELET-TRANSFORMED IMAGES 
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Fig. 8. — Raw images of line emission (O continuum, Mg XI, Si XIII, S XV, Ar XVII, Ca XIX, and Fe xxv) in Cas A and corresponding 
wavelet-transformed images for five different scales. The white scale bar is 1' fa 1 pc in length. The color bar is set so blue is the minimum, 
and red is the maximum. 




Fig. 9. — Raw images of line emission (Fe L, Mg XI, Si XIII, S XV, Ar XVII, Ca XIX and Fe xxv) in Kepler and corresponding wavelet- 
transformed images for five different scales. The white scale bar is 4:l"Ki 1 pc in length. 



Si XUl 



S KV 




Fig. 10. — Raw images of line emission (Fc L, Si XIII, and S xv) in Tycho and corresponding wavelet-transformed images for five different 
scales. The white scale bar is 87"~ 1 pc in length. 



sxv 




Fig. 11. — Raw images of line emission (Si XIII and S xv) in G15.9+0.3 and corresponding wavelet-transformed images for five different 
scales. The white scale bar is 50"~ 1 pc in length. 
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Fig. 12. — Raw images of line emission (Mg XI, Si XIII, S xv) in Gil. 2— 0.3 and corresponding wavelet-transformed images for five different 
scales. The white scale bar is 41"Ri 1 pc in length. 
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Si xui 




Fig. 13. — Raw images of line emission (Mg XI, Si XIII, S xv) in Kes 73 and corresponding wavelet-transformed images for five different 
scales. The white scale bar is 52"Ri 2 pc in length. 
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Kes73 




Fig. 16. — Twenty-three substructures in Gil. 2— 0.3 (left), thirty substructures in Kes 73 (middle), and thirty-one substructures (right) 
identified with WTA. We extracted X-ray spectra at these locations, and the fit results are given in Tables 4-6. 
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Fig. 17. — Example Chandra X-ray spectra, models, and residuals from Gil. 2— 0.3 (left; from circle Q), Kes 73 (middle; from circle F), 
and RCW 103 (right; from circle R). 



A.2. DETAILED SPECTRAL MODELING 

Since three of the nine SNRs in our WTA sample (Gil. 2— 0.3, Kes 73, and RCW 103) do not have pubhshed Chandra 
X-ray spectra and modehng of their ejecta, we provide here a more detailed analysis and presentation of the spectra 
from these targets. 

We identified 23 substructures with WTA in Gil. 2— 0.3. We extracted Chandra spectra for the seven ACIS obser- 
vations of Gil. 2— 0.3 from regions centered on the 23 substructures. Since the goal of this analysis is to describe the 
physical state of the ejecta, here we chose to improve the statistics by increasing the radii of the regions where we 
extracted spectra to 30 pixels = 14.76"(the cyan circles A-W in Figure [T51 left). Background spectra were produced 
from a region sa50" from Gil. 2— 0.3 and subtracted from the source spectra. Spectra were modeled like the analyses 
above, with an absorbed plane-parallel shocked plasma model, and data from the seven observations were fit simulta- 
neously to improve statistics. We let the abundances of magnesium, silicon, sulfur, and iron vary freely with the other 
elements frozen to solar values in the fits. Example spectra and models from one region (circle Q in Figure 1161 left) 
are given in Figure [T7] (top) . 

Table 4 lists the parameters of the best-fit models. The mean absorbing column density is high, (Nh) = 2.4x10^^ 
cm~^, and attenuates the soft X-rays below about 1 keV. The electron temperature kT varies from 0.60 keV up to 
1.38 keV and has a mean value of (kT) — 0.80±0.18 keV. The Northern (regions C and D) and the Western (regions 
I and J) regions seem to be hotter than other parts of the remnant, and the area just south of the pulsar (regions 
U, V, and L) appear to have elevated temperatures as well. The abundances of magnesium, silicon, and sulfur are 
supersolar, consistent with a shocked ejecta origin of the X-ray emission lines. The only exception is region W, which 
is near solar metallicity, suggesting that the X-ray emission is from shocked ISM. The ionization timescale n^t, the 
product of ambient electron density and the time since the plasma was shock heated, spans an order of magnitude, 
with Uf-t = 9.47 X 10^° — 1.04 x 10^^ s cm~'^. The mean ionization timescale from our fits is {riet} = 4.85 x 10^^ s cm~^. 

In Kes 73, we identified 30 substructures with WTA. We extracted the Chandra spectra at these locations with 
radii of 30 pixels = 14.76"(cyan circles A-d in Figure [TOl middle). Background spectra were produced from a region 
«50" from Kes 73 and subtracted from the source spectra. Spectra were modeled as above, with a single plasma in 
non-equilibrium ionization. We let the abundances of magnesium, silicon, sulfur, and iron vary freely with the other 
elements frozen at solar values during the fits. An example spectrum and model from one region (circle F in Figure [TCI 
middle) are given in Figure [171 (middle). 

Table 5 lists the parameters of the best-fit models. As in Gil. 2— 0.3, the mean absorbing column density is high, 
(Nh) = 2.7x10^^ cm~^, and attenuates the soft X-rays below about 1 keV. The Nh varies across the remnant, from 
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TABLE 4 
Spectral Results for Gil. 2-0. 3 
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2.1-3.4x10^^ cni^^, and appears to be elevated in the northern portions of Kes 73 (e.g., regions C and D). The electron 
temperature kT ranges from 0.63 keV up to 1.22 keV and has a mean value of {kT) = 0.84±0.49 keV. Generally, the 
abundances of magnesium, silicon, and sulfur are supersolar, with a few exceptions: seven regions (R, S, T, a, b, c, 
and d) are consistent with solar metallicity within their errors. This result suggests that the southern and Eastern 
sections of Kes 73 are dominated by X-ray emission from the shocked ISM. The ionization timescale Uet spans nearly 
two orders of magnitude, with riet = 7.75 x 10^" — 5.24 x 10^^ s cm""^. The mean ionization timescale from our fits is 
{net) = 7.85 x 10" s cm'^. 

In RCW 103, we identified 31 substructures with WTA. We extracted the Chandra spectra at these locations with 
radii of 30 pixels — 14.76"(cyan circles A-e in Figure [TBI right). Region O was larger, with a radius of 60 pixels since its 
Omax was large, w60 pixels. Background spectra were produced from a rectangular region with wl50xl5" sides south 
of RCW 103 (chosen to avoid chip gaps) and subtracted from the source spectra. Spectra were modeled as above but 
excluding energies above 3 keV because of a dominant non-thermal component. We let the abundances of magnesium, 
silicon, and iron vary freely with the other elements frozen at their solar values during the fits. An example spectrum 
and model from one region (circle F in Figure [TO]) are given in Figure [T71 (bottom). 

Table 6 lists the parameters of the best-fit models. These fits were successful statistically, except for one region, 
circle O in Figure I16[ which was poorly fit by one NEI plasma. The addition of a second NEI plasma improved that 
fit statistically, with a Ax^ = 58, for 126 degrees of freedom. The mean absorbing column density of our RCW 103 
spectral models is moderate, {Nh) = 5.4x10^^ cm~^, so the soft X-rays are relatively unattenuated and the broad 
Fe L emission is strong. The electron temperature kT is fairly constant across the remnant, with all models giving 
kT K, 0.52-0.67 keV except one, region L, at 0.33 keV. The mean electron temperature is {kT) — 0.56±0.14 keV. 
The iron abundances have roughly solar values, suggesting the Fe L emission comes from the shock-heated ISM. The 
magnesium and silicon have supersolar abundances, with only a couple outliers, so the Mg XI and Si xiii lines likely 
have an ejecta origin. The ionization timescales of RCW 103 are generally greater than those of Gil. 2— 0.3 and Kes 
73, consistent with an older age for this source. The mean ionization timescale from our fits is {nj) — 9.25 x 10^^ s 
cm^^, indicating the plasma in RCW 103 is approaching collisional ionization equilibrium. 
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TABLE 5 
Spectral Results for Kes 73 
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122/110 


a 


2, 


.3 


0.88±0.08 





.87±0.19 


1, 


,07ib0, 


.19 


1, 


.00±0.21 


2, 


.14cll 


123/106 


b 


2, 


.1 


0.89ib0.11 


0, 


.98ib0.24 


1, 


.18ib0, 


.22 


1, 


.12ib0.26 


2, 


.lOell 


98/87 


c 


2, 


.3 


0.84ib0.10 





.92ib0.28 


1, 


.28ib0, 


.32 


1, 


.09ib0.30 


3, 


.28ell 


80/79 


d 


2, 


.2 


1.22ib0.16 


1, 


.13ib0.33 


1, 


.llibO, 


.24 


1, 


.37ib0.31 


1, 


.09ell 


110/95 



X-ray Morphological Properties of SNRs 



TABLE 6 
Spectral Results for ROW 103 



Region 



(xlO- 



Nh 

22 ^„ 



kT 
(keV) 



Fe/Fe0 Mg/Mg. 



Si/Sio 



n^t x^/d-o.f. 



(s cm 



A 
B 
C 
D 
E 
F 
G 
H 
I 
J 
K 
L 
M 
N 
O 
P 

Q 

R 
S 
T 
U 
V 

w 

X 
Y 
Z 
a 
b 
c 
d 
e 



0.6 
0.8 
0.6 
0.5 
0.6 
0.7 
0.7 
0.7 
0.6 
0.8 
0.7 
1.0 
0.5 
0.6 
0.5 
0.5 
0.4 
0.5 
0.5 
0.5 
0.3 
0.3 
0.4 
0.4 
0.4 
0.4 
0.5 
0.4 
0.5 
0.5 
0.4 



0.52±0.03 
0.53±0.03 
0.57±0.01 
0.55±0.01 
0.52±0.01 
0.52±0.03 
0.61±0.03 
0.67±0.03 
0.53±0.02 
0.60±0.02 
0.57±0.02 
0.33±0.06 
0.55±0.02 
0.59±0.02 
0.56±0.01 
0.58±0.02 
0.55±0.01 
0.57±0.01 
0.53±0.02 
0.60±0.01 
0.63±0.03 
0.58±0.01 
0.55±0.01 
0.59±0.01 
0.54±0.02 
0.53±0.01 
0.57±0.01 
0.59±0.02 
0.57±0.02 
0.66±0.01 
0.54±0.01 



1.35±0.23 
0.98±0.15 
1.05±0.13 
0.97±0.10 
0.88±0.08 
1.65±0.20 
1.09±0.17 
1.52±0.29 
1.28±0.16 
1.09±0.19 
1.11±0.14 
1.58±0.44 
0.92±0.26 
0.91±0.12 
0.81±0.03 
0.62±0.06 
0.67±0.05 
0.89±0.08 
1.24±0.14 
0.83±0.08 
0.56±0.24 
0.73±0.07 
0.76±0.06 
1.00±0.09 
0.91±0.09 
1.02±0.10 
0.85±0.07 
1.20±0.13 
0.55±0.08 
1.04±0.10 
l.llitO.ll 



1.75±0.26 
1.19±0.17 
1.47±0.16 
1.42±0.13 
1.41±0.11 
1.84±0.22 
1.70±0.23 
1.86±0.28 
1.51±0.18 
1.84±0.27 
1.40±0.16 
1.51±0.31 
1.72±0.37 
1.31±0.15 
1.12±0.04 
1.02±0.10 
1.08±0.08 
1.15±0.10 
1.43±0.15 
1.06±0.10 
1.47±0.46 
1.03±0.09 
1.00±0.08 
1.29±0.11 
1.14±0.11 
1.34±0.12 
1.18±0.09 
1.50±0.15 
1.23±0.14 
1.75±0.15 
1.75±0.16 



2.19±0.37 
1.69±0.25 
1.68±0.21 
1.38±0.17 
1.56±0.16 
1.84±0.26 
1.12±0.21 
1.40±0.27 
1.82±0.25 
1.66±0.28 
1.12±0.17 
2.94±1.20 
1.40±0.43 
1.41±0.19 
1.28±0.06 
1.03±0.12 
1.13±0.11 
1.08±0.11 
1.82±0.22 
1.05±0.13 
1.83±0.68 
1.11±0.13 
0.98±0.11 
1.42±0.15 
1.21±0.16 
1.18±0.15 
1.45±0.13 
1.41±0.18 
1.21±0.18 
1.52±0.16 
1.78±0.21 



5.61ell 
6.43cll 
7.49cll 
1.88cl2 
1.88cl2 
2.70cll 
5.17cll 
2.67cll 
4.73cll 
7.24cll 
5.76ell 
6.83ell 
2.04el2 
6.73ell 
1.07el2 
1.54el2 
1.57el2 
7.70ell 
7.57cll 
1.09cl2 
1.71cl2 
9.08cll 
1.03cl2 
4.19cll 
4.69cll 
6.58cll 
7.95cll 
4.68ell 
2.06el2 
5.34ell 
9.05ell 



60/63 
74/88 
116/79 
98/88 
116/91 
126/92 
122/86 
126/91 
102/88 
105/88 
86/76 
72/68 
122/82 
114/100 
348/126 
112/108 
145/98 
148/110 
111/99 
126/91 
53/69 
147/105 
131/107 
153/101 
88/96 
133/100 
127/97 
127/84 
80/77 
115/93 
86/87 



